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Abstract. 

When a ligand that is bound to an integral membrane receptor is pulled, 
the membrane and the underlying cytoskeleton can deform before either the 
membrane delaminates from the cytoskeleton or the ligand detaches from the 
receptor. If the membrane delaminates from the cytoskeleton, it may be further 
extruded and form a membrane tether. We develop a phenomcnological model for 
this processes by assuming that deformations obey Hooke's law up to a critical 
force at which the cell membrane locally detaches from the cytoskeleton and a 
membrane tether forms. We compute the probability of tether formation and show 
that they can be extruded only within an intermediate range of force loading rates 
and pulling velocities. The mean tether length that arises at the moment of ligand 
detachment is computed as are the force loading rates and pulling velocities that 
yield the longest tethers. 

E-mail: tomchou@ucla.edu 

1. INTRODUCTION 

Adhesion between cells plays an important role in a number of biological processes 
involving cell motility and cell-cell communication. Cell-cell adhesion is mediated by 
integral membrane proteins on the surfaces of interacting cells. Cadherins, which 
mediate binding between cells of the same type within a tissue, bind to themselves, 
while integrins, which mediate binding between different cell types bind to Inter- 
Cellular Adhesion Molecules (ICAMs) or Vascular Cell Adhesion Molecules (VCAMs) 
[I] . Understanding the physics of this adhesive interaction requires an understanding 
of both the protein-protein bond as well as the cell's mechanical response when these 
bonds become stressed. If forces act to pull the two cells apart, either of the cells' 
cytoskeletons and plasma membranes may deform. Under certain conditions, the lipid 
membranes can delaminate from the underlying cytoskeleton and be pulled into long 
tethers. At any time during this process, the bonds holding the two cells together may 
also break, arresting tether extraction. 

A specific biological process in which bond dissociation and membrane 
deformation must be considered simultaneously is leukocyte extravasation, which is 
part of the process by which leukocytes are recruited to inflamed or infected tissue. 
Endothelial cells that make up blood vessels preferentially express cellular adhesion 
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molecules, including selectins, near wounded tissue. Leukocytes circulating in the 
blood can then bind to the endothelial cells via their own cell surface proteins. Bonds 
between the leukocytes and the endothelial tissue are transiently made and broken as 
a shear flow in the blood vessel pushes the leukocyte, rolling it across the endothelial 
layer [2]. The rolling leukocytes contain microvilli that are enriched in adhesion 
molecules that preferentially attach to endothelial cells. During rolling, these microvilli 
tethers can extend under the hydrodynamic shear force of blood flow in the vessel [3] . 
At the same time, forces imposed on the endothelial membrane via the adhesion 
molecule can cause the endothelial membrane to form a tether [3]. Tether formation 
and extension of microvilli can decrease the force that the adhesive bond feels from 
the shear blood flow. 

Both the physics of cell membrane deformation and the mechanics of ligand- 
receptor bonds have been studied extensively. Micropipette, atomic force microscope 
(AFM), optical trap, and magnetic bead techniques have been used to pull membrane 
tethers from cells and probe the properties of the cell membrane [5]. The diameter of 
the extracted tether depends on the membrane surface tension and bending rigidity 
and, if the tether is being extended at a constant velocity, the membrane viscosity [HE]. 
Therefore, these quantities can be inferred from tether pulling or poking experiments 
[8]. Theoretically, Euler-Lagrange methods have been used to compute equilibrium 
tether shapes and the force-extension curve for a tether pulled quasi-statically from a 
lipid vesicle [9j[10]. These theoretical models of pure lipid bilayers show only a ~ 10% 
overshoot, or barrier, in the force-extension curve before a tether is extracted from 
an asymptotically flat membrane [9j liojfl . some experiments show a significant force 
barrier to tether formation [Til Q21 H3] . These large force barriers to membrane tether 
formation arise in living cells and is attributed to membrane adhesion to the underlying 
actin cytoskeleton [UJ[T4]. When tethers are pulled from giant artificial vesicles, the 
size of the force barrier can increase only when the area on which the pulling force is 
exerted is increased [12]. For smaller vesicles, area and volume constraints may also 
influence the tether force-extension relationship jl5] . These may arise from nonlocal 
terms in the functional describing the lipid membrane energetics. For example, area- 
difference elasticity can give rise to a restoring forces that continually increase as 
tether length increases [16] [17] . Such nonlocal effects will be important only when the 
tether comprises an appreciable fraction of the total membrane area. Force curves that 
do not saturate at long tether extensions can also arise when part of the membrane 
reservoir adheres to a substrate [13j . The relative importance of nonequilibrium forces 
arising from the viscosity of both the membrane lipids and surrounding solution have 
also been estimated [18]. 

During tether pulling experiments, and in their corresponding theoretical models, 
the molecular bond attaching the pulling device to the membrane is assumed to 
always remain intact. However, a typical ligand-receptor bond used to connect the 
pulling device to the lipid membrane (and possibly the underlying cytoskeleton) can 
rupture upon pulling. Although the details of a bond's energy landscape can be 
probed using dynamic force spectroscopy [HI [20] [21] [22], one can usually assume 
that bond rupturing is dominated by a single activation barrier that can be lowered 
by an externally applied pulling force. Most AFM studies of molecular strength 
are performed on proteins that have been isolated from cells. However, recent 

f The barrier represents the energy that must be overcome before a lipid tether can be drawn out. 
The ~ 10% energy barrier arises from the Helfrich free-energy-minimizing geometry of an incipient 
lipid membrane tether in the absence any cytoskeletal attachments |91 [TOl (cf. dashed curve in Fig. [2] 
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studies have probed bond strengths between proteins still embedded in a live cell 
membrane [23] HI] • While using live cells has the advantage that the post-translational 
modifications of membrane proteins are preserved, the mechanical deformation of the 
cell's cytoskeleton and membrane must also be taken into account. To model this 
system, a viscoelastic Kelvin model was used to fit experimental measurements of the 
force-extension relationship to determine effective cellular adhesion [24] . 

In this paper, we develop a dynamic model that incorporates phenomenological 
theories of membrane and cytoskeleton deformation, tether extraction, and the kinetics 
of ligand-receptor detachment. In contrast to an equilibrium model determining tether 
extraction and detachment from an adhered vesicle [15] , we find relationships that 
define when tether extraction is likely, and the typical length of the tether pulled 
before the ligand-receptor bond ruptures. In the next section, we motivate a simple 
mechanical model using phenomenological forms for the force-extension relationship of 
the membrane. Given the large bending energies of lipid bilayer membranes, and the 
relatively strong attachment of membranes to cytoskeleton, we will neglect thermal 
fluctuations of the membrane, but implicitly include thermally-driven ligand-receptor 
bond dissociation. Dynamical equations are written for two commonly employed 
experimental protocols, a linear force loading rate, and a constant bond pulling speed. 
In the Results and Discussion, we compute the probability of tether formation and 
plot universal curves that delineate regimes where tethers are likely to form. Mean 
tether lengths are also plotted for both pulling protocols. 



2. Mathematical Model 



Consider the system depicted in Figured] A ligand is bound to an integral membrane 
protein, which may also be directly associated with the cytoskeleton. The ligand may 
be attached to a pulling device via a cantilever spring and is pulled with either a fixed 
speed, or a force that increases linearly in time. 

As the ligand is pulled, the cytoskeleton first deforms, and eventually can detach 
from the membrane. At this point, the lipid membrane may flow into a tether. At 
any point during the cytoskeletal or membrane deformation, before or after membrane- 
cytoskeleton delamination and tether formation, the ligand affixed to the pulling device 
may detach from the membrane receptor protein. 



2.1. Membrane Mechanics 



We first consider the response of the membrane-cytoskeleton system to an externally 
applied pulling force F p (t). The rate at which the receptor-ligand complex moves will 
be described by 



dx(t) 
dt 



dU(x) 



dx 



F P (t) 



C (t) 



(1) 



where £ is a mobility that is inversely related to the viscosity of the membrane lipid 
[T8] . In general, although £ depends on the configuration of the system defined 
predominantly by x(t), we neglect the details of this dependence and assume it to 
be constant. 

The term U(x) represents the energetic cost associated with deforming the cell 
membrane and underlying cytoskeleton when the receptor is displaced a distance x 
normal to the flat membrane. This phenomenological energy can be derived from 
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Figure 1. (a) A cell membrane can be pulled by a device such as a cantilever. 
The device can be moved at fixed velocity V, or, at a specified force F p (t). Under 
specified force conditions, the transduction of force through the device is assumed 
instantaneous such that F p (t) is always acting across the ligand-receptor bond. 
Before membrane delamination from the underlying cytoskeleton, the pulling 
results in small membrane and cytoskeletal deformations, (b) After membrane 
delamination, the pulling device can extrude long lipid tethers without deforming 
the cytoskeleton. The total system length z is the sum of the displacement x of 
the membrane protein from its initial position and £, the increase in the length of 
the cantilever spring from its unstretched length. 



a detailed consideration of the membrane-cytoskeletal mechanics. For simplicity, we 
will assume the membrane mechanics are governed by a Helfrich free-energy |25j that 
includes a lipid bilayer bending rigidity and an effective thermally-derived entropic 
membrane tension. We will assume that the membrane reservoir is large enough such 
that an extruded tether negligibly depletes the reservoir. Hence, global contributions 
to the membrane energetics, such as area-difference elasticity, can be neglected. 

Experiments in which tethers were pulled from live cells found a significant force 
barrier to tether formation [TTJ [14] . While a smaller force barrier can also arise in 
pure lipid membranes [9] , we will assume that the plasma membrane is attached to an 
underlying cytoskeleton (with anchoring molecules) , which we model as a linear elastic 
material provided the deformation is small. The receptor that binds the ligand that is 
attached to the pulling device can also be a transmembrane receptor that is directly 
attached to the cytoskeleton. As the membrane is initially pulled, the cytoskeleton 
will elastically deform as a Hookean spring. The receptor or anchoring molecules 
will break at a deformation xq, and the lipid membrane will be drawn into a tether. 
This occurs at a critical delamination force F c . Thus, for displacements x < Xq, 
where the Hookean approximation for the membrane-cytoskeleton assembly is valid, 
the membrane-cytoskeleton carries an effective spring constant F c /xq, where F c is the 
critical delamination force. Experimentally, the cytoskeleton typically detaches from 
the membrane where the filament-free tether forms [26J [27] ; therefore we can assume 
a linear force extension relationship of the form 
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At a displacement x = x (h the membrane delaminates from the cytoskeleton and 
a lipid membrane tether forms. Under the infinite reservoir and simple Helfrich 
free-energy assumption, the tether can then elongate indefinitely under a constant 
force dll/dx = Fq that is intrinsic to the lipid tether and is determined by the 
membrane bending rigidity k and entropic surface tension a [9j [10] : Fo — 2tt\ / 2k<t. 
The phenomenological membrane force-displacement relationship for both attached 
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Figure 2. Wc assume that when the receptor is initially associated with the 
cytoskeleton, d x U(x) increases linearly until some maximum force F c in the barrier 
to tether formation is reached, after which the tether extends with constant force 
-Fo (solid curve) . The qualitative features of this force-extension curve is observed 
in numerous systems [7] 1111 1121 1131 I14| . For reference, we show d x U(x) as a 
function of x when the only contributions to U (x) come from membrane surface 
tension and bending rigidity (dashed curve). This curve was calculated using the 
method described in [9] assuming a membrane bending rigidity of 20k bT and a 
surface tension of 0.0076 dynes/cm. 



(solid curve) and free (dashed curve) membranes are shown in Figure [2l The barrier 
F c to tether formation is larger for receptors or membranes that are attached to the 
underlying cytoskeleton, than for the free lipid membrane case. 

In order to close the equations of our basic model, we must specify F p (t). 
Henceforth, we will consider two cases typically realized in experiments: a linearly 
increasing (in time) pulling force, and a fixed pulling speed. 



2.2. Linear Force Ramp 



For a force linearly increasing in time, F p (t) = Tt, where T is the rate with which the 
force increases. Eq. [T] becomes 



dx(t) 
dt 



dU{x) 
dx 



Tt 



(3) 



Upon defining the time t at which the membrane delaminates from the cytoskeleton 
provided the ligand- receptor bond has not yet ruptured by x(t — t a ) = x , we have 
(x/x )F c for t < t Ql and §^ = F for t > t . Thus, Eq. |3]is solved by 



dx 
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and 

x(t>t )= x -tF (t-to) 

where the delamination to is determined from the solution to 
^o-^d-e-^o )=1 . 



(5) 



(6) 



2.3. Constant Pulling Speed 

In the case of constant pulling speed, we must include the dynamics of the device 
deformation £(t). Since the velocity of the pulling device is fixed, we note that the 
total displacement z{t) — x(t) + £(t) obeys 



dz dx d£ 
dt ~ ~dl + dt ~ 



(7) 



This equation holds only when the ligand is attached to the membrane-bound receptor. 
Let us assume that the pulling device has an internal response that is modeled by a 
simple spring so that the force F p (t) that the spring exerts on the ligand is 



F p (t) = K£(t), 



(8) 



where K is the spring constant of the pulling device. Since the pulling force is 
proportional to £(t), it will ultimately depend on the pulling rate V and the physical 
properties of the pulling device (represented by an elastic cantilever in Figure 1) and 
cell membrane through Eq. [7J Upon integrating Eq. [7] and using the initial conditions 
x(t = 0) = lit = 0) = 0, we find 



s(t) = Vt-£(t). 



(9) 



Substituting Eqs. O and [5] into Eq. [TJ and expressing the dynamics in terms of the 
device deformation, we find a closed equation for £(t): 



d£ 

dt= V + t 
This equation is solved by 

£(t < t ) = 



dU(x) 



d.r 



K£(t) 



x=Vt-i(t) 



(10) 



VKx _0 ^ _ e -(,{K+F c /x )t^ 



C(F C + Kxo) 2 



for t < to, and 



£(t > to) 



F c Vt 
~F C + Kx 



_ V + £F _ , K{t _ to) 
~ UK { ' 

+^(t )e-« K( *- to) 



(11) 



(12) 
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for t > t (when dU/dx = F ). Here, the time to at which tether formation occurs 
is found by evaluating Eq. [HI at time t = to, x(t ) = x — Vt — t(to), yielding an 
implicit equation for t : 



v* FcVt ° 
Vt - x 



F c + Kxo 



VKx l _q _ e -(L{K+F o / X0 )t y 



(13) 



i{F c + Kxo) 2 



After evaluating to numerically, £(t) is found in terms of the K, V , £, F Cl Fq, and xq, 
and the membrane displacement can be found using Eq. [9] 

2.4- Ligand- Receptor Dissociation 

The dynamics described above for the membrane and pulling device deformations 
assume that the pulling device remains attached to the membrane through an unbroken 
ligand-receptor bond. Since all external forces are transduced through the ligand- 
receptor bond, the pulling force F p (i) on the membrane (cf. Eq. [1} vanishes once 
the ligand-receptor bond ruptures. However, the probability of ligand-receptor bond 
dissociation itself depends on the applied force F p (t). We can model the breaking 
of the ligand-receptor bond by a Poisson process and define a ligand-receptor bond 
survival probability Q(t) that obeys 

^ = -k r (t)Q(t), (14) 

where k T (t) is the force-dependent rupture (or dissociation) rate of the ligand from 
the receptor. We assume that /c r (i) takes a simple Arrenhius form |28) : 

fc r (i) = fc e Fp(t)d/fcBT , (15) 

where d is the length of the ligand-receptor bond and ksT is the thermal energy. The 
solution to Eq. [M] is explicitly 



Q(t) = exp 



fc / e ^(*'WfeBT dt / 



(16) 



More complex models of dynamics bond rupturing can be derived [29l [30]. Here, 
for simplicity, molecular details such as the thermally-induced bond-breaking attempt 
frequency and the intrinsic free energy of the unstressed ligand-receptor bond are 
subsumed in the effective rate parameter fco- 

Since w(t) = —dQ(t)/dt — k r (t)Q(t) is the bond rupture time distribution, 
the mean membrane displacement at the time of ligand-receptor rupture (the mean 
maximum displacement) is given by 

/>oo 

w(t)x(t)dt= / k r (t)Q(t)x(t)dt. (17) 
./o 

In our subsequent analysis, we will combine bond rupturing statistics with 
membrane tether dynamics and explore, as a function of the physical parameters, 
the probability of tether formation, and the length of pulled tethers should they form. 
To be concrete, we will use typical parameter values (listed in Table I) found from the 
relevant literature to guide our analysis. 
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parameter 


range of values 


reference 


d 


0.8-1.0 nm 


31 (streptavidcn-HABA) 


£ 


~ l//m/(pNs) 






~ 1 — 4/zm 


SH 


Po 


3-380pN 


calculated^ 




100-380pN 


[27] (red blood cells) 




l-100pN 


[32] (epithelial cells) 


K 


8-llpN//im 




k a 


10- b -10/s 


33 


V 


~ 3 /im/s 




r 


lOpN/s 


n 



3. Results 

Here, we compute the dynamics of tether formation and ligand-receptor bond 
rupturing under both linear force loading and constant pulling velocity protocols. 



3.1. Linear Force Ramp 



When F p (t) = Tt, the ligand-receptor survival probability defined by Eq. [16] is 
explicitly 



Q(i) = exp 



k a k B T 
Yd 



{e 



Ttd/k B T 



1) 



(18) 



The bond survival probability, evaluated at the time of tether formation to (found 
numerically from Eq. Q(to) = Pt, determines the likelihood that a tether is 
extracted, and is plotted in Fig. G^a) as a function of force loading rate T. Note 
that Pt first increases with the force loading rate T, before decreasing again at very 
high loading rates. Large critical delamination forces F c increase the probability 
that ligand-receptor bonds detach before membrane-cytoskeleton delamination occurs. 
Membrane-cytoskeleton combinations that have weaker delamination forces F c yield 
a larger range of force loading rates that lead to tether formation. Moreover, since 
the pulling force is specified, Q(t) is independent of the free tether restoring force Po 
(cf. Eq.lTS"]). The only dependence is on the delamination force F c which sets the 
delamination to (cf- Eq.[S]) in the expression Pt = Q{to)- pnce the tether is formed, 
the free tether restoring force Fo is irrelevant. The values T + and T_ are defined by 
Px(r±) = 1/2 and define the window of loading rates within which tether formation 
is likely. 

For quantitative evaluation of T±, and their dependences on the other system 
parameters (P ma x, d, ksT, xo, and £), it is convenient to define dimensionless 
parameters according to 

_ gr _ k x d £,F C 

7=T2 — i a =77r^F' /c=j , (19) 

kpXo t,k B i k x 

and find parameter regimes within which Pt > 1/2. Upon using Eq. 1181 the phase 

boundaries for tether formation are determined from the implicit solution to 

1 



exp 



1 



(20) 



t F = 2V2tt^k^ where re is the membrane bending rigidity and a is the effective membrane surface 
tension[9]. We assumed cr = 3 - 1200pN//im[34] and re = 10 - 20k B T [35] 





Figure 3. (a) Tether formation probability Pr = Q(*o) as a function of 
force ramp rate F Note that results for the constant force ramp protocol are 
independent of the free tether restoring force Pq- r± denote the linear loading 
rates at which the tether formation probability Pp = 1/2. Applied loading 
rates T_ < T < r_|_ are likely to lead to tether formation. For F c = 20pN, 
r_ PS 7.2pN/s and S3 435pN/s. Other parameters used were: ko = 0.01/s, 
d = lnm, xo = lfim and £ = l/xm/(pNs). (b) Dimensionless parameter regimes in 
which Pp > 1/2. The regions of parameter space below each curve are associated 
with Pr > 1/2, where tether formation is likely. The smaller the dimensionless 
bond dissociation rate a = koxod/ {(.kgT), the wider the range of dimensionless 
loading rates 7 = ^F/(kgXo) leading to tether formation. The maximum and 
minimum pulling rates <y™ in = 7— (/ c = 0) and 7™ ax = 7+(/ c = 0) are indicated 
for the a = 0.006 curve, while the maximal dimensionless delamination force 
jmax j s snown f or a = 0.004. (c) The minimum and maximum force ramps, and 
the maxima delamination force as functions of a. For 7 < 7™ 111 , 7 > 7™ ax , or 
f c > /™ ax , tethers always have less than 50% chance of forming. 



where r = k to is determined from the solution to the dimensionless form of Eq. [5] 
^-^(l-e- /cT0 ) = l. (21) 

Jc Jc 

Figure [3jb) shows curves in dimensionless parameter space (/ c and 7) below which 
_Pt > 1/2. Asymptotic analysis of the condition Pt = 1/2 shows that for sufficiently 
large 7, tether formation will always be suppressed. Additionally, as the dimensionless 
ligand-receptor dissociation rate a increases, the regime for tether formation shrinks. 

Conditions for tether formation can be further refined by computing universal 
parameter curves that define regimes for which Pt can never be greater than 1/2. As 
a function of the intrinsic ligand-receptor dissociation rate, there is a band of pulling 
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rates outside of which Pj is always less than one half, even when tether extraction 
is barrierless (/ c = 0). Fig. [3fc) shows 7™ 111 and 7? lax , the minimum and maximum 
dimensionless ramp rates at which P? = 1/2 when f c = 0. These delimiting loading 



rates are roots of Eq. For small a, the lower root 7" 



2/ ln^ 2 + 0(a) 



Also plotted in Fig. |3jc) is /,f ax (a), the maximum membrane-cytoskeleton 
delamination force that can give rise to Pt > 1/2, for any ramp rate. One is unlikely 
to pull tethers from membranes that require more force than / ( ? lax to delaminate from 
the cytoskeleton. Moreover, there is a critical a max « 0.22444 above which = 1/2 
cannot be reached, independent of / c or 7. 

Finally, we plot in Fig. [4] the expected dimensionless maximum tether length 
(X*) = (x*)/xo found from Eq. [TTl as a function of the dimensionless force loading 
rate 7. Not only does the mean tether length decrease with increasing critical 
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Figure 4. (a) Mean tether length (X*) as a function of load rate 7 for different 
critical delamination forces f c . (b) {X*) at fixed f c and various /o. Even though 
the cytoskeleton-free membrane restoring force does not affect tether formation 
probability Px, it does influence the length of tether possible. 



delamination force f c , as shown in Fig. |4|a), but the value of 7 at which (X*) is 
maximized decreases as / c is raised. Since tether extrusion is a competition between 
the rate of climb against the restoring force F c x/xq and ligand- receptor dissociation, 
as f c is increased, higher ramp rates are required to cross the delamination barrier 
faster, relative to the bond rupturing. Although the free membrane restoring force 
/o = £,Fo/(koXo) does not come into play in the tether formation probability, it does 
influence the mean length of lipid tether that is extruded before the ligand detaches 
from the membrane-bound receptor. Note that in the constant loading rate protocol, 
the force is specified and all results are independent of the pulling device rigidity K. 



3.2. Constant Pulling Speed 



Now consider a constant pulling speed protocol. The force felt by the membrane in 
this ensemble will depend on the will depend on the pulling device rigidity K . The 
bond survival probability is computed from 



Q(i) = exp 



-k 



,Ki(t')d/k B T dt > 



(22) 
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where t(t) is given by Eqs. llll and [T2"l Initially, while the membrane and cytoskeleton 
are attached to each other, and constant speed pulling is applied, the ligand-receptor 
bond survival probability Q(t) first decreases rapidly. After delamination, the force 
on the ligand-receptor are fixed, arising only from Fq and the viscous drag The 
subsequent decay of Q(t) arises from a slower, single exponential. 

Fig. EJa) shows the corresponding Pj as a function of pulling speed V, and as in 
the force ramp case, reveals an optimal pulling speed that maximizes the likelihood 
of pulling a tether. In the V — > limit x(t) sa (from Eqs. l9l and fT2|) . and we expect 
Pt — > because when the ligand-receptor bond detaches, the membrane has not been 
sufficiently deformed. In the fast pulling limit, the detachment rate increases quickly 
and ligands detach at very short times t such that tether formation cannot occur. 
Upon defining Pt(V±) = 1/2, pulling speeds V— < V < V+ are likely to result in 
tethers. 
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Figure 5. (a) The probability Pt that tether formation occurs before the ligand 
detaches from the receptor plotted as a function of pulling speed V. Analogous 
to the force ramp protocol, there is a band of pulling speeds V— < V < V+ 
within which the tether formation probability is greatest. Parameters used were: 
k = 0.01/s, d = lnm, £ = l^m/(pN s), F c = 20pN, K = 10pN//»a. (b) Phase 
diagram in the f c = £F C / (koxo) and v = V/(koxo) parameter space for various 
a = koxgd/ (kgT) and fixed pulling device stiffness /i = ^K/kg = 1000. (c) Phase 
diagram for fixed a = 0.002 and various pulling device stiffnesses fi. 



To explore how V± depends on other system parameters, we employ the same 
dimensionless parameters defined in Eqs. 1191 along with a dimensionless pulling speed 
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and pulling device stiffness, 

V K£ 

v = and /i = — — . (23) 

k x k 

Figure [S^b) shows, for /i = 1000 and various values of a, the boundaries below 
which tether formation is likely. Fig. [5jc) shows the phase boundaries for fixed 
a = 0.002 and various pulling device stiffnesses /i. Note that softer pulling devices 
suppress tether formation at low speeds since the forces are not immediately felt by 
the membrane, allowing the ligand more time to detach. However, softer pulling 
devices greatly enhance tether formation at large pull speeds because the accelerated 
delamination more than compensates for the drag-mediated acceleration of ligand- 
receptor dissociation. 
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Figure 6. (a) Universal curves for the critical values v™ ln and i;™ ax as a function 
of dimensionless unstressed ligand-receptor dissociation rate a. Curves for three 
difference values of device stiffness fi are shown. For a <C 1 and /i — > oo, the 
minimum pull speed is asymptotically u™ m ~ l/(ln2 — a). In this limit, the 
critical a max beyond which u™ ln and u™ ax merge (and tether formation is unlikely, 
even for f c = 0) occurs at a max « 0.24, slightly larger than the a max in the 
constant load rate protocol (not shown), (b) The maximal delamination force 
jmax a f unc tion of a for two different stiffnesses fi. 

Note that analogous to the linear force ramp protocol, for particular a and 
H, there is a maximum and minimum pulling speed (u™ m and u™ ax ) beyond which 
tether formation is unlikely, even when the delamination force vanishes (Fig. UK a)). 
Conversely, there is a maximum delamination force /™ ax above which no pulling speed 
will result in likely tether formation (Fig. EJb)). 

Finally, consider the mean receptor displacement at the moment of ligand 
detachment, (X*) = (x*}/x , found from Eq. [17] with X{t) = vt — £(t)/xo, and 
the appropriate Q(t) and k T (t). The same arguments that explain the non-monotonic 
behavior of Pt as a function of V apply here, and the mean dimensionless receptor 
displacement at the moment of ligand detachment, (X*), is a non-monotonic function 
of V (Fig. [7]). When the pulling velocity is small, the receptor moves slowly and the 
term x(t) in Eq. [T7] will be small, rendering integral in Eq. [T7] small. On the other 
hand, when V is sufficiently large, the viscosity £ -1 allows a large force to be reached, 
accelerating ligand-receptor bond rupturing. Thus, Q(t) quickly decreases and the 
mean receptor displacement when the ligand detaches, (x*), will be small. 




Figure 7. (a) The mean system length at ligand detachment, {X*), is a non- 
monotonic function of the dimensionless pulling velocity, v, and increases with 
decreasing f c . The qualitative dependence of {X") on fo is similar to that shown 
in Fig. [4] (b) The mean maximum tether length as a function of v for various 
dimensionless pulling device stiffnesses /j. 



Both constant force load rate and constant pulling speed protocols give 
qualitatively similar tether extraction probabilities and maximum delamination forces 
/™ ax . This is not surprising since in the two protocols are physically equivalent in 
both the infinitely stiff and infinitely soft pulling device limits, prior to delamination. 
In the large \x oc K limit, a constant pulling speed forces the ligand to have the 
trajectory x(t) = Vt. Because we assume a linear force-extension relationship before 
delamination, this protocol is equivalent to a high constant force loading rate of 
r « KV (or 7 s» jiv). In the extremely soft pulling device limit, the elastic pulling 
device absorbs most of the extension and the force at which it acts on the ligand also 
increases linearly in time: T sa F c V/xq (or 7 sa vf c ). 

However, we do find a qualitative difference in the mean tether length extracted, 
due to difference in the post-delamination forces between the two protocols. As 
functions of load rate and pulling speed, the maximum mean tether lengths attainable 
via linear force ramp are typically less than half of those achieved through constant 
pulling speed, all else being equal. This feature can be understood by considering 
how the receptor displacement, x(t) and the ligand dissociation rate, k T (t) depend on 
time in each case. When the pulling speed is constant, after the tether forms, x(t) 
increases linearly in time, and k r (t) is constant. When we apply a force ramp to the 
system, x(t) increases quadratically in time once tether formation occurs. This would 
seem to imply longer tethers under the force ramp protocol; however, in this case, the 
dissociation rate k r (t) also increases exponentially in time. Thus, ligand detachment is 
much faster in the force ramp case, resulting in shorter observed mean tether lengths. 

4. Summary and Conclusions 

We modeled membrane-cytoskeleton delamination in series with a ligand-receptor 
bond and a deformablc pulling device and determined the parameter regimes within 
which lipid tether extrusion is likely. Results from our model can be directly used to 
propose and analyse experiments in which cell or lipid vesicle membranes are pulled by 
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a breakable bond. For example, in [TT], tethers are pulled from endothelial cells when 
large force barriers are overcome, but detachment of the pulling device from the tether 
is not considered. Performing such experiments with breakable ligand-receptor binds 
would provide the necessary data with which to test our predictions on the likelihood 
of tether formation and on the differences between fixed load rate and pulling speed 
protocols. 

For both linear force ramp and constant pulling speed protocols, we find a wide 
window of ramp rates and pulling speeds that likely lead to tether extraction. However, 
we also find critical values of a dimensionless membrane-cytoskeleton delamination 
force, and a dimensionless spontaneous ligand-receptor dissociation rate beyond which 
tether formation is unlikely, regardless of all other parameters. We assumed in all of 
our analysis that the tether force-extension curve can be derived from local interactions 
with a Helfrich free-energy model. Finite-size membrane reservoirs and nonlocal 
energies such as area- difference elasticity would give rise to increasing forces as the 
tether is extended, thereby increasing the probability of ligand-receptor dissociation, 
and decreasing expected tether lengths (X*). 

Both linear force ramp and constant pulling speed protocols yield intermediate 
tether formation regimes, with a specific pulling speed v and specific linear ramp rate 
7 that maximizes the mean tether length (X*) in the respective protocol. However, 
they present different tether dynamics after delamination leading to different expected 
tether lengths (x*). Using both protocols, and our results, it may be possible to 
characterize membrane-cytoskeleton properties, provided sufficent information about 
the ligand-receptor binding energy and pulling device response are known. In general, 
such inverse problems are very ill-posed, but restricting the force-extension relationship 
to simple forms as we have done, one may be able to use the onset of tether formation 
as a way to estimate force parameters (such as Fq,F c and ieo)- 

While we have framed our analysis in terms of AFM experiments in which the 
strength of a ligand-receptor bond is probed while the receptor is in the membrane of 
a live cell, our basic model is relevant to leukocyte rolling as well. In leukocyte rolling, 
a bond between protein on a leukocyte microvilli and a protein in the membrane of 
an endothelial cell becomes stressed. Because the microvilli act like Hookean springs 
when pulled [3J, and simultaneous extension of microvilli and tether extraction from 
the endothelial cells has been observed [4], our analysis is directly applicable to this 
system, with the cantilever replaced with a microvilli. 

Finally, we note that we have treated the ligand-receptor bond rupturing as a 
stochastic Poisson process, while the deformation of membrane and cytoskeleton was 
considered deterministic. This approximation is good as long as xq ^> d. However, 
if the experiment is repeated, each region of membrane may have highly variable 
attachments to the cytoskeleton. In this distribution of delamination forces 

F c should be considered. Another source of stochasticity may arise when multiple 
adhesion points are being pulled, possibly leading to multiple tethers [35] ■ If the entire 
system is treated as a single, effective tether, the force-extension of this super-tether 
will rely on the statistics of how many individual tethers are still attached during the 
dynamics. 
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